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Abstract. FEARLESS (Fluid mEchanics with Adaptively Refined Large Eddy Simu- 
lations) is a numerical scheme for modelling subgrid-scale turbulence in cosmological 
adaptive mesh refinement simulations. In this contribution, the main features of this 
tool will be outlined. We discuss the application of this method to cosmological simu- 
lations of the large-scale structure. The simulations show that the production of turbu- 
lence has a different redshift dependence in the intra-cluster medium and the warm-hot 
intergalactic medium, caused by the distinct stirring mechanisms (mergers and shock 
interactions) acting in them. Some properties of the non-thermal pressure support in 
the two baryon phases are also described. 



1. Introduction 

In the framework of the physics of the cosmological large-scale structure, turbulent gas 
flows are an important link between the thermal and merger history of galaxy clusters, 
on the one side, and the non-thermal phenomena (cosmic ray acceleration, amplifica- 
tion of magnetic fields) in the intra-cluster medium (ICM) on the other. A central role 
in injecting turbulence in the cosmic flow is played by shocks, which contribute both to 
energy dissipation and gas stirring ( Miniati et al.,.2000: Vazza et al..,2()09.). and by hv- 
drody namical instabilities triggered by mergers (e.g.. lHeinz et al. llapichino et al.1 
2008|) . 



The evolution of turbulence in the ICM, as a result of the cluster merger his- 
tory, has been explored in many w o rks over the last decade using hydrodynamical 



simula tions (Ricker & Sarazin 200 F; "bolag et all 120051 : llapichino & Niemeved 120081 



Vazza et al.,t2011. : ,Paul et al.,>2011.) . In this contribution we wan t to highlight the result s 



of a different approach to the study of turbulence, adopted by l lapichino et al.l ( 2011 ). 



Instead of focusing on a single cluster, that work followed the evolution of turbulence in 
a large cosmological box (with side length of 100 Mpc h~^). This setup allows to study 
not only the injection of turbulence in the ICM, but also in the less dense warm-hot in- 
tergalactic medium (WHIM). Furthermor e, the simulation co de includes a subgrid scale 



(SGS) model for unresolved turbulence (ISchmidt et al.ll2006) . coupled to the adaptive 
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mesh refinement (AMR). The resulting tool, called FEARLESS (Fluid mEchanics with 
Adaptively Refined Large Eddy Simulations), will be briefly outlined in the next sec- 
tion. 



2. Numerical methods 



FEARLESS consists of the combination of AMR with a SGS model for the unre- 
solved kinetic energy. Details of this numerical tool have been presented elsewhere 
dSchmidt et al.ll2006l : lMaier et al.l2009l : llapichino et al.l20Tl '): in essence, the discretiza- 
tion ont o a grid of the equations of fluid dynamics is equivalent to applying a. filter for- 
malism dGerm ano 1992) to them. Consequently, additional terms appear in the equa- 
tions, which take into account the dynamics at unresolved length scales. For example, 
the filtered momentum equation of a viscous, compressible, self-gravitating fluid, be- 
comes 



dt 



dr: 



dr: 



dr: 



dr: 



(1) 



where p is the baryon density, v, are the velocity components, p is the pressure, gi the 
gravitational acceleration and cr'.. the viscous stress tensor. Given a variable /, with / 
we indicate the application of the filter operator to it (see ISchmidt et al.ll200^ for the a 
more rigorous treatment). 

The last term on the right-hand side of equation ([T]) contains the turbulent stress 
tensor f(v;,vy), which accounts for the interaction between the resolved flow and the 
SGS scales. This term can be expressed in analo gy wit h the viscous stress tensor by 
means of the so-called eddy viscosi ty closure (cf. Pope 2000l). a lthough an improved 
closure has been recently adopted by Schmidt & FederrathI ( 201 1 ). 

The turbulent stress tensor enters also the definition of the specific filtered kinetic 
energy, as a contribution from unresolved scales: 



^kin = ^V/V,- -I- ^T(v,-,Vy)/<p) 



(2) 



It is thus natural (lGermanolll992 ) to interpret the trace of t(v,-, Vj)/{p) as the square of 
the SGS turbulence velocity q. This leads to the definition of the SGS turbulence energy 
et. 

1 , 1 

(3) 



The SGS turbulence energy is governed by an equation of the following form: 
-{p)e, + ^Vj{p)er - D + X + F - {p){A + e) , 



dt 



(4) 



The quantities on the right-hand side of equation dUl determine the evolution of Ct and 
are the turbulent diffusion term D, the turbulent production term E, the pressure di- 
latation term A and the viscous dissipation term e. Their closures represent the SGS 
model. 

The method is coupled with AMR in a way that consisten tly accounts for cut-off 
length scales varying in time and space (see iMaier et al.ll2009l for a detailed descrip- 
tion). The resulting tool is particularly suitable in the study of turbulence in strongly 
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Figure 1 . Temporal evolution of the mass-weighted averages of the internal (eint) 
and turbulence SGS (et) energies for the ICM and WHIM baryon phases. From 
llapichino etaTI (1201 ll) . 



clumped media. Furthermore, the cell-wise computation of et makes the analysis and 
visualization of turbulence easy and flexible. 

FEARLESS has been im plemented on the gr id-ba sed, AMR hybr i d (jV-b ody plus 
hydrodynamical) code Enzo dO'Shea et all 12005 ). In llapichino et al.l ( 201 ih . we ran 
adiabatic simulations in a cosmological box with a side of 100 Mpc h~^, resolved on 
a root grid of 128^ cells and 128^ A/^-body particles. Four additional AMR levels have 
been allowed, with a refinement criterion based on overdensity, leading to an effective 
spatial resolution of 48.8 kpc h~^. 



3. Turbulence and non-thermal pressure support 



The definition of the baryon phases in llapichino et al.l ( 2011 ) is based on a threshold in 



temperature T and baryon overdensity 6. The gas is defined as belonging to the WHIM 
if r > 10^ K and ^ < 10^; if 5 > 10^ it belongs to the ICM. The former phase is found 
in filaments and cluster outskirts, and the latter in the denser, collapsed structures. In 
Fig. [T] the time evolution of the average internal and SGS energy for these two phases 
is shown. We notice a different redshift evolution for et. in the ICM, it shows a peak 
at z between 1.0 and 0.65, followed by a mild decrease. In the WHIM phase, there is a 

steady increase to z = . 

In ll apichino et al.l ( 2011 ) we provide an interpretation of these features in terms 



of different stirring mechanisms acting in the two baryon phases: the evolution of tur- 
bulence in the ICM follows closely the cluster n ierging history, ther efore et peaks ap- 
proximately during the major merger epoch (cf. iGiocoli et al.ll2007l) an d its d ecline is 
halted by the subsequent minor mergers. Recently, Hallman & Jeltema (| 201 lb studied 
the evolution of turbulence in clusters, and found that the fraction of clusters with large 
turbulence in the core evolves in time with a trend very similar to the ICM in Fig. [T] The 
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evolution of turbulence in the WHIM is governed by the gas accreted on filaments and 
cluster outskirts, and closely resembles t he evolution of th e kinetic energy flux through 

external shocks dMiniati et al. 2000: Ski lhnan et al.ll2008h . 

Another interesting problem which was explored in llapichino et al.l (1201 lb con- 
cerns the dynamical pressure support of the cosmic gas. We notice that, starting from 
the definition of ei, one can identify the non-thermal pressure caused by unresolved, 
SGS velocity fluctuations with pt - 2/3 pet. This term has been included i n an anal- 
ysis o f the dynamical supp ort against g ravitational contraction of the gas (IZhu et al. 



20101) . We refer to lapichino et all (1201 ih for a more thorough derivation of the support 



equations; here it is sufficient to say that the analysis is based on the Laplacians of the 
thermal and dynamical pressure (the latter one including both resolved and SGS terms). 

It is found that the turbulent support is stronger in the WHIM gas at baryon over- 
densities I < 6 < 100, and less relevant for the ICM. A fairly large fraction of the 
WHIM and ICM gas has a large vorticity (28.7 and 52.3 per cent in mass at z = 0, 
respectively), but this is usually associated with an equally large thermal pressure sup- 
port, and only in a small volume fraction (of the order of 10 percent) the non-thermal 
pressure is dynamically significant. 

This result is apparently in contradiction with the idea that the cluster outskirts, 
consisting mostly of gas which is newly accreted in the potential well, could have sig- 
nificant departures from hydrostatic equilibrium. O ur work shows th at this is the case 
only in localized regions, but not globally (see also IValdaminil 12011 ). More resolved 
cosmological simulations, focused on single clusters, will be used to investigate further 
this point. 
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